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1. INTRODUCTION 

The study of voltage behavior in power grids has become a significant concern for the operators and 
planners of these systems. Several widespread incidents around the world have been associated with voltage 
instabilities. This mode of uncertainty is not yet well-controlled, compared to the method of angular instability 
(dynamic and transient). Indeed, the mechanism causing voltage instability seems to be one of the most critical 
issues to be clarified. Even today, there is no widely available and widely accepted theory for analyzing voltage 
stability. The problems linked to this type of instability then constitute, in several countries, a central axis of 
research. Power systems are vulnerable to voltage instability due to unpredictable power demand. The strategy 
requires that the system remain in a secure operating state under abnormal conditions [1]. The electric power 
systems are stressed caused by an increased and fluctuating loading. This leads to the overloading of 
transmission lines, voltage instability, and could contribute to voltage collapse blackouts [2]. In fact, blackouts 
in power networks are mainly the result of a miscellaneous process that involves several factors [3]. 
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The instability of tension plays a major role in the occurrence of these extraordinary events. Some important 
identation examples of these incidents are reported in [4-8]. 

The blackout in Denmark and Sweden reported on september 2003 [4]. The system disturbance that 
affected Algeria power network in February 3, 2003 [5]. Another major blackout was reported on mid August 
2003 that affected the northeastern and midwestern of United States, and the Canadian province of Ontario. 
The high reactive output power, due to the contact with trees, caused the initial disturbances [6]. The disturbance 
that affected the electrical power network of Italy system in September 28, 2003 [7]. The Moscow’s system 
disturbance of 25 May 2005, the Brazil’s system disturbance of November 10, 2009, the Tokyo system 
disturbance of March 11, 2011, and the India’s system disturbance of July 31, 2012 [8] are few examples. 

New monitoring techniques make it possible to improve the control of transits on existing 
transmission lines and increase network transport capacities (loadability). These technological solutions are 
based on power electronics are generally called FACTS, an acronym of "Flexible Alternative Current 
Transmission Systems". The FACTS was born in the 1970s, and certain technologies are currently installed on 
the electrical network. Thus, the SVC (static var compensator) and their improvements, the STATCOMs 
(known as static synchronous compensator) are the most mature technologies. The goal of their development 
is to allow better voltage control in normal operation and in case of faults. The optimal positioning of FACTS 
devices in a network is a combinatorial problem. To this day, there is no analytical method that can solve this 
kind of problem and give the overall optimum. Therefore, tools such as the metaheuristic algorithm have been 
used. Evolutionary algorithms (EA), such as the particle swarm, genetic algorithm, and differential evolution, 
have been commonly used in the power system. These algorithms were originally used for optimizing 
single-objective problems. In addition, evolutionary methods provide more flexibility in handling 
the complexities in objective functions and constraints. Poorly known just a few years ago, optimal positioning 
and sizing of FACTS devices is considered to improve electrical network loadability. Enhancing network 
loadability is crucial in modern electrical grids. FACTS devices are used to reduce short circuit currents and to 
enhance both transient and steady-state stability. Several approaches have been proposed for the optimal 
placement of the FACTS*s devices using different conventional and metaheuristic optimization methods. 

The authors in [10] proposed two approaches to identify the placement of the FACTS controller of 
the suitable nodes to ensure voltage stability. The first technique is based on the topological structure of power 
networks, while the seond technique is based on the conventional power flow [10]. In [11], a multi-objective 
allocation problem of FACTS devices is solved by two-hybrid approaches; 1) non-dominated sorting particle 
swarm optimizer (NSPSO) combined with the fuzzy logic, 2) non dominated sorting genetic algorithm 
II (referred as NSGA-ITI) combined with the fuzzy logic. In this work, power loss, L index voltage stability, and 
voltage deviation are optimized simultaneously. The goal of [12] is to find the optimal location and setting 
parameters of SVC and thyristor controlled series capacitor (TCSC) devices, using PSO to mitigate small signal 
oscillations in multi-machine power systems.The authors in [13] present the application of particle swarm 
optimization to determine the optimal position of the FACTS devices in the most cost-effective manner. 

The optimal placement of STATCOMS is suggested in [14].Through the simultaneous application of 
PSO and CPF in order to improve the voltage profile to efficiently minimize power losses and to maximize 
system loadability. The continuation power flow (CPF) method allows calculating the static voltage stability 
margin when only the STATCOMSs size is considered. The PSO and the general algebraic modeling system 
(GAMS) were also adopted to obtain location and rating of D-STATCOM for voltage stability margin 
enhancement [15]. In the literature, CPF has been used to solve the problem of the singularity near the stability 
limit of the conventional power flow algorithms. The authors in [16] investigate the voltage stability assessment 
using static VAR compensation (SVC) and thyristor controlled series capacitor (TCSC) controllers. The effect 
of SVC and TCSC sizing in the loading margin was also evaluated [17]. 

The aim of this paper is to determine the optimal placement of the most popular FACTS devices, 
namely the SVC. The PSO objective function which contains a term that needs to be maximized is the loading 
factor. The optimization procedure was executed under one and three SVC constraints. In the first step, we 
consider only the reference voltage constraint and in the second step, three constraints are taken into account, 
the reference voltage, the reactance and the SVC reactive power in the purpose to maximize the load level and 
minimize the total losses power system and flattening the voltage of the buses. To the best of our knowledge, 
this is the first study that addresses the problem of finding an optimal location of the SVC s device under those 
optimization constraints. The remainder of this paper is organized as follows: section 2 introduces the voltage 
stability analysis. Particle swarm optimization is presented in the section 3. The problem is formulated in 
section 4. The simulation results are analyzed and discussed in section 5. Finally, our contributions and 
conclusions are given in section 6. 
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2. VOLTAGE STABILITY ANALYSIS 
2.1. Bifurcation parameter 

Bifurcation theory is concerned with the study of topological type changes in the trajectories solution 
of a non-linear dynamic system when the control parameters vary. It focuses mainly on the evolution 
of equilibrium periodic orbits (which constitute asymptotic states), their multiplicity and their stability 
properties [16-18]. Parameter estimation is a complementary function of the state estimator necessary to ensure 
its robustness and reliability. These parameters are used in several advanced functions of network control, such 
as for example: contingency analysis, monitoring of safety limits, optimal power flow, etc. In this theory, 
the system equations depend on these parameters, as (1). 


fQ.) = 0 (1) 


The parameter (A) is used to reproduce the load changes that drive the power systems to voltage collapse. 
The load power is modified as: 


Pri = Prio(1 + AK;;) (2) 


Qiri = Qrio(1 + ÀK) (3) 


frepresents the power flow equations, y is dependent variables, (A) is the loading parameter, Prion 
and Qio represent the active and reactive powers at the basic operating point respectively. K,; interpret 
the changing of active and reactive power at bus i as (A) varied. In a typical bifurcation diagrams, the voltages 
are represented as functions of the system loadability measurement (A). This representation is called P-V 
or nose curves, as shown in (2), and they are widlely used in the continuation power flow (CPF) analysis. 


2.2. Continuation power flow CPF 

Power systems P-V curves and the maximum loading parameter are calculated using extensively CPF 
techniques. Thus, an iterative process is used. This process in a discrete context comprises 2 stages, 
the prediction (predictor step) and the update (corrector step). The idea is to obtain a better estimate for each 
sample of time t. Figure 1 illustrates the iterative process; from a known initial position A, also called solution A, 
a tangent predictor is used to estimate a new solution for a specified pattern of load increase (called solution B). 
While the system load assumed to be fixed, a corrector step (update iteration) 1s used to determined 
the exact solution (or solution C) by using a conventional power flow analysis. This process is iterated until 
the desired P—V curve is obtained. In order to optimize the calculated loading parameter, we have developed 
an optimization algorithm by combining the CPF and PSO, which is described in the following section. 


Predictor 


7 I 
. B Corrector 


Exact 
Solution 


Bus voltage 












Critical voltage 


Loading Parameter (A) 


Figure 1. An illustration of the CPF technique 


3. PARTICLE SWARM OPTIMIZATION 

Swarm intelligence (SI) refers to a distributed computing paradigm based on collective intelligence 
that emerges from the cooperation of multiple autonomous agents. The SI applies concepts based on the fodder 
of ant colonies and the grouping of social animals such as birds using agent-based calculations. The three major 
paradigms of Swarm Intelligence are: “Ant Colony Optimization’, “Particle Swarm Optimization”, and robotic 
swarms (“Swarm Robototics“ or SR) [19]. The optimization by a swarm of particles is based on a set of 
individuals originally arranged randomly and homogeneously, which we will henceforth call particles, which 
move in the hyperspace of research and constitute, each a potential solution. Each particle has a memory 
concerning its best-visited solution as well as the ability to communicate with the particles constituting its 
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surroundings. From this information, the particle will follow a tendency made, on the one hand, of its will to 
return to its optimal solution, and on the other hand, of its mimicry compared to the solutions found in its 
vicinity. From local and empirical optimums, the set of particles will normally converge towards the optimal 
global solution of the problem treated. The PSO is a substitute to genetic algorithms and ant colonies for 
the optimization of non-linear functions [20]. 

Binary string coding is a modification of the PSO algorithm to solve problems with solution elements 
with binary values. The meaning of the speed variable has been changed to indicate the probability 
of the corresponding solution element with a value of O or 1. The speed is updated in the same way as for 
the classic PSO. No coefficient of inertia is used. The performance of each particle (i.e., its proximity to 
the global optimum) is measured by means of an objective function, which varies with the optimization 
problem [20]. As with other swarm algorithms, choosing the right parameters is important for efficient PSO 
execution. Considerable work has been done to select a combination of values that work well across a wide 
range of issues. Let’s consider y and v as the coordinates and velocity of the particle, respectively. The particles 
move to take into account their best position and best vicinity (panurgic displacement), this will repeated for 
each iteration. The position of each particle in the next step is then evaluated as the sum of its current position 
and the speed is represented as pBest. 

The index of the best particle is represented as gBest; then the swarm can be carried out by solving 
the following update equations: 


Va+1 = W; ' Va + Q1 ' rand(-) * (pBest — vq) + p2 ' rand(-) * (gBest — Xa) (3) 
Xat1 = Xa t Vas (4) 


whered is the iterations index, y,is the current particle position at the d iteration. vg is 
the particle velocity at the d‘ iteration. W; represents the inertia weight factor, @, and p, are acceleration 
constant. rand(-)is uniformly distributed random number in the interval [0,1], ite%, q, represents 
the maximum number of iterations while iter is the number of the iterations until the current stage. Proper 
selection of the inertia weight W; provides a balance between the global and the local exploration. The W; 
is set according to (5). Additionally, Figure 2 illustrates the PSO particles behavior. 

Wi max—Wi min x Iter (5) 


Itermax 


Wi = Wimax — 


At first glance, it seems that several parameters have been considered for the PSO algorithm. However, 
several of those parameters can be fixed in advance; others, on the contrary, can only be defined empirically. 
This is the case with the swarm population size. There is no rule to determine this parameter; doing many tests 
allows you to gain the experience necessary to understand this parameter. We must also consider 
the initialization of the swarm; it is generally done randomly according to a uniform distribution on [0,1]. In 
the above procedure, the particle velocity is limited by the maximum value vmax. Two other important 
parameters are the confidence coefficients previously named @, and ø, (also called acceleration constant). 
They make it possible to balance the tendencies of particles to follow their conservation instinct. Similarly, an 
important parameter to take into account is the coefficient of inertia (Wi) in the formula seen before. It defines 
the exploration capacity of each particle in order to improve the convergence of the method. Setting this 
parameter amounts to finding a compromise between a global exploration (Wi>1) and a local exploration 
(Wi<1). It represents the adventurous instinct of the particle [21-23]. Thus, a suitable selection of the inertia 
weight Wi using (5) allows achieving a better exploration of the search space. 


Ug 


gBest 





EE 
Xa a 
© pBest 


Figure 2. The PSO particles behaviour 
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4. PROBLEM FORMULATION 
4.1. Objective function 

Throughout the previous section, we have seen the theoretical aspects of swarm intelligence, and more 
specifically, PSO. The usual application of PSO remains optimization. From the moment when the semantics 
of the problem to be solved can be expressed in the form of a function to be optimized, the PSO applies. 
In this paper, an optimization process is used is to obtain an efficient utilization of the existing power network 
with SVC optimal location. In this scope, the cost function that need to be maximized can be written as: 


Maximize (A) (6) 


4.2. SVC modeling 
The time constant regulator is assumed to be known using the model described in Figure 3. 
The absolute reactance bsvc is well defined. The following differential equation can be used [24]: 


; POD_y)\_ 
Dege = (Kr(Vref +Vs v) Vsyc) (7) 


Tr 


The total reactive power introduced at the SVC device is expressing as follows [24-25]: 
q = bsyc V’ (8) 


where V is bus voltage magnitude, VP?” is the signal output of the power oscillation damper, V,.. fis reference 
voltage, K, is regulator gain, T,„is regulator time constant, bmax is maximum susceptance, bmin 1s Minimum 
susceptance and bsyc is total reactance. The optimization process maximizes (à) while satisfying the following 
inequality constraint. 





vV Hn 


Figure 3. The structure of the SVC 


4.3. Equality constraints 
The typical load flow equations can be given by: 


Poi = Ppi = V; D Vi (Gi; cos(6; = 6;) F Bij sin(d; = 6;)] = 0 yl = 1, ng NB (9) 
Qgi — Qpi a V; 2 V| Gi; sin(d; — 6;) F Bij cos(6; = 6;)| = () yl = 1, siny NB (10) 


where Pęand Qg are respectively the real and reactive power generator. Pp and Qpare respectively 
the real and reactive power load. G;jand B;; are respectively the transfer conductance and susceptance between 
bus i and bus j while NB represents the number of buses of the electric network. 


4.4. Inequality constraints 
An optimization problem with constraints is a combinatorial problem where we want to maximize an 
objective function. We, therefore, find ourselves in an optimization context where it is not only a question 
of satisfying a set of constraints but also of achieving a so-called optimal solution. These constraints used in 
this work represent the system operating and they are represented by: 
- Electricity generation constraints: the volatage Vc and reactive power outputs Qg are limited by their lower 
bounds V2", Q@™" and upper bounds Ve~, Q7°*, respectively, where NG is the number of generators. 
and they can be described as follows: 


Ve eV, TH dag NG (11) 
OO. T= le NG (12) 
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- Transformer constraints: The transformer ratios are discrete values although they are here modeled as 
continuous variables. Transformer tap T settings, which are based on transformer rating and system voltage, 
where NT is the number of transformers, are bounded as follows: 


Oe a. i= Ayo NE (13) 


- SVC constraints: These include the constraints of reference voltage V,¢,, the total reactance bsyc and 
the reactive power Qsyc, where NSVC is the number of SVCs as follows: 


Vere <Vrepi < Veer, i= 1,...,NSVC (14) 
pre eh br, T= 1, ime SVC (15) 
e020, d= wig NSVC (16) 


The proposed PSO-CPE algorithm is described in the flowchart of Figure 4. 





Initialize PSO parameters 


Determine and store pBest and gBest 


Figure 4. Flowchart of the proposed PSO-CPF 


5. CASE STUDIES 

In this formulation, the objective function remains quadratic, but it is unweighted. Rather, standard 
deviations are used as variables that add to the measures. In addition, binary variables are integrated in order 
to detect erroneous parameters. A PSO-CPF algorithm based on the improvement of the objective function is 
proposed to help solve this mixed variable problem. For our simulation results, we will use the PSO-CPF 
algorithm that we developed. To improve the voltage stability of the electric network, an optimal location of 
the SVC is considered. Different scenarios will be tested for the IEEE-30 node network. Our simulation steps 
will go through 2 cases, in the first the constraint is taken only the voltage of the SVC, and for the second case, 
the voltage of the FACTS device thus considered its reactance and its reactive power in the objective of 
maximizing the loading parameter taken as the objective function. 


5.1. Simulation result for IEEE 30-bus system 


First, the formulation is tested on the IEEE 30-bus network. Errors ranging from 5% to 100% of 
the initial parameter values are introduced and correctly identified. Even in cases where the susceptance and 
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transformer ratio errors are adjacent, the parameters can be detected and estimated with accuracy. The load 
cases (or Operating points) are represented by the real and reactive powers necessary for each bus in 
the network. Some bars linked to a generator generate power. Other buses are connected to subnets and require 
a certain outgoing power load. The power flows in the branches adjust according to the generation and 
the demand between the bus [26]. Figure 5 illustrates the network with the data considered.The susceptance 
error varies from 1% to 100% of its initial value. The measurement of reactive power flow in the branches is 
where susceptance has the most impact. The redundancy of this measurement, as well as its amplitude, varies 
according to the lines. Table 1 presents the cases tested and PSO parameters. The conclusions which can be 
drawn from them. 


Table 1. The IEEE-30 Bus system PSO 




















Parameter name PSO variables 
Population size 30 
Weighing factor Wi 0.9-0.4 
Constant ol 1.4 
Constant 2 1.4 
The number of iterations 100 
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Figure 5. The IEEE-30 bus system 


5.2. Case 1: SVC voltage reference only 

In this scenario, the voltage reference of the SVC was taken into account. The optimal solutions 
obtained from the PSO and their related results are provided in the Table 2. Based on the voltage profiles of 
Figure 6, at the critical buses (buses 9, 10, 11, 12, 14, 15, 16, 17, 18, 19, and 20) the voltages results obtained 
have increased with SVC. A weak destabilization occurs at buses 5, 7, 8, 26, 29, and 30 because the SVC has 
not been installed on the weakest bus. The active power losses profiles of the system (without and with 
the presence of the SVC device) are shown in Figure 7, near the collapse point the increase in losses is higher 
in the case 1 of the SVC device. Figure 8 shows the reactive power losses profiles of the system for the two 
cases without and with the presence of the SVC. We can see the increase in losses near the collapse point is 
lower in the observed of the case 1 of the SVC. Figure 9 illustrates the improvement of the objective function 
according to the iterations. It is interesting to note in this figure that the convergence characteristics of the 
proposed PSO algorithm. 


Table 2. The best results for case 1 
The best location of SVC The best à (pu) The best Vref (pu) 
12.0000 -3.2519 0.9500 
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Figure 6. 30-bus system voltage profile without and Figure 7. 30-bus system active power losses profile 
with the presence of the SVC device (case 1) without and with the presence of the SVC device (case 1) 
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Figure 8. 30-Bus system reactive power losses Figure 9. The convergence characteristics of 
profile without and with the presence of the proposed PSO algorthim applied to 
the SVC device (case 1) 30-Bus system (case 1) 


5.3. Case 2: all SVC constraints 

In this case, the whole set of problem constraints were taken into consideration. The optimal solutions 
obtained from the PSO are summarized in Table 3. The voltage profile of our simulation is shown 
in Figure 10. On the other hand, the Figures 11 and 12 show the active and the reactive power losses profile, 
respectively. The best position for this case was founded is bus 15. For this scenario the corresponding size 
of SVC is 62.105 MVAR, -16.98 MVAR (inductive) to +45.125 MVAR (capacitive). It can be seen that 
an smart and optimal SVC placement could reduces both the real and the reactive power losses. This represents 
a considerable gain for voltage stability. 

At the collapse point the SVC device provides a better voltage profile as illustrated in Figure 10. 
This is because the SVC is installed at the weakest bus, which is the bus 15. From the active and reactive power 
losses profiles of the system (without and with the presence of the SVC device), we noticed the increase in 
losses near the collapse point is lower in the case 2 of SVC.The convergence characteristic of the proposed 
PSO for case 2 is depicted in Figures 13. Figure 14 illustrates the impact of installing the SVC devices. 
This confirms the previous observations. Therefore, it is clear from this figure that total losses are minimized 
after the SVC installation on bus 15. 


Table 3. 30-Bus system results (case 2) 


The best location of SVC The best à (pu) The best Vref (pu) The best b (pu) 
bmin bmax 


15 3.1823 0.9500 -0.1882 0.5000 
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Figure 13. 30-bus system 
convergence haracteristics of 
the proposed PSO (case 2) 


Figure 12. 30-bus system reactive power losses 
profile without and with the presence of 
the SVC device (case 2) 





14 

12 

10 
8 
6 
4 : ; 

B Without SVC 
2 
E With SVC Case 1 
g 
" ~ E With SVC Case 2 
D 4 
C a a rr 
ee se SC we 
of gh 
E & Ra gè E 
é a a +2 k 
A oes g ms 7 „a 
5 Oo añ fe 
AI Pi oF = a qe 
ms i A se 40 = 
ab a a a 
ao si a 
wt 
40 


Figure 14. The impact of the SVC 


6. CONCLUSION 
In this paper, a new formulation of the parameter estimation function for an electrical network has 


been presented. The resulting model is a combinatorial optimization problem, and a PSO-CPF combination 
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algorithm has been proposed to estimate the optimal placement of SVC, and therefore to enhance the power 
system voltage stability. The formulation was tested on the IEEE 30 bar network with the Gaussian noise of 
the measurements. However, since this new formulation has been developed in transport networks where noise 
is high, this detection of small errors is not an objective. The new formulation is therefore competitive with 
the other methods of the literature review in addition to being independent of the weights of the database. Other 
observations and potential improvements have emerged from the various tests carried out; the reactive powers 
of the branches need to be evaluated to avoid false detection. This work can be further extended to address 
the problem of the optimal location of FACTS devices and wind generators to enhance voltage stability. 
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